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Abstract 

The relativistic equation of motion of a particle with spin is derived. The 
effects of the spin are compared to the effects of the self force. For low 
energies the spin effects are shown to be two orders of magnitude larger than 
the self force. For higher energies is is shown that the self forces are bigger, 
and the overall magnitude of the effects are compared. 



1 Introduction 

With reports of laser intensities of 10 22 W cm -2 , and expectations of an increase by 
at least an order of magnitude, pQ the classical equation of motion of electrons in 
electromagnetic fields is gaining renewed interest. When considering the motion of 
electrons in such extreme fields, these particles quickly become free and attain 
energies that vastly exceed their rest energies. The correspondence principle, with 
such extremely high principle quantum numbers, clearly indicates the validity of the 
classical approximation, where not only radiation reaction effects become 
important, [2] but the effect of the interaction of the intrinsic magnetic moment of 
the particle and the external electromagnetic field becomes large. [3] 

An unsolved problem relates to the comparison of the intrinsic spin effects to 
the self self forces. Specifically, we need to understand at what energies spin effects 
dominate and if, and when, self forces predominate. A simple calculation in the low 
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velocity limit indicates that the spin effects are important. Let us consider the ratio, 
R a , self force to the spin force in the low velocity limit, 



2e 2 d/3 



(1) 



V(/x B 

and use a = eE/m, where the dot indicates the time derivative, e is the charge, a is 
the acceleration, \x is the magnetic dipole moment, and c = 1. For an 
electromagnetic wave of the form Ecos(kx — out) this yields 



Ro = a (2) 

where a is the fine structure constant. This result shows that the spin force is two 
orders of magnitude larger than the radiation reaction force. Of course this result 
holds in the low velocity limit, and indicates that the spin force should be examined 
more carefully. If one naively (and incorrectly) replaces da/dt by da/dr in the above, 
than the quotient in (j2J) is multiplied by v°, which becomes large for high energies. 
This indicates that at high energies the self forces may become larger than the spin 
forces, and a more careful analysis in the following will show that this is the case. 

The point, and an objective of this paper, is to find the relativistic 
generalization of this result. This will be achieved by first developing the classical 
equation of motion with spin, ignoring radiation reaction effects. Then the equation 
of motion with self forces will be examined, ignoring spin, and the results will be 
compared. In this approximation we ignore the radiation reaction effects of the spin 
part, which is a small part of a small part. 

An early detailed work on this subject is given by Bjabha and Corben.[l] 
This subject was also treated in the book of Barut,[5] who returned to the problem 
of radiation reaction years later. [6] In the 1960s, a classical equation of motion with 
a spin interaction was derived by starting with the quantum mechanical equation, 
and obtaining the classical limit by disregarding the commutators. j7j A year later, 
following an attempt by Pauli, an asymptotic solution to the Dirac equation was 
developed using the WKB method for a particle with an anomalous moment, [8] and 
an anomalous moment was used to explore classical electrodynamics resulting from 
a nonlinear Dirac equation. [9] A comprehensive review of particle-laser interaction 
in high intensities has been published recently. |10j 
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This article derives the equation of motion from an action principle, but also 
derives the field equations that follow from the same action. The generalized BMT 
equation is also derived. A few definitions will help clarify the discussion below. 
The spin tensor is defined in terms of the spin vector S a according to 

Sfj,u — ^fiuapS 01 !) 13 (3) 

where e^ va ^ is the totally antisymmetric tensor. Equivalently 

S, = VVe^. (4) 
Other useful definitions are the magnetic moment tensor, 

//* = —S a (5) 
m 

and the electromagnetic dual "F^ u = \F a pe al3 ^ v . 

Finally, the traditional equation of motion of a charged particle in an 
electromagnetic field is given by 

f = eF^ (6) 

where p M = rraA 

2 Equation of motion with spin 

We know that the force on a magnetic dipole in a magnetic field is given by 

/ = V(/i.B). (7) 

This well known result is derived by integrating the Biot-Savart law around a small 
current loop and defining the magnetic moment as a line integral in the usual way. 
It is not a valid definition for an elementary particle, since we know that the 
magnetic moment cannot be explained by a small loop current. Therefore we 
introduce the interaction directly into the action. The relativistic generalization of 
the scalar fi ■ B is 



(8) 



which implies the generalization of (JTJ) is = (e/2m)(F fll ,S^ u )'' i and therefore the 
naive relativistic generalization of ([H]) is dp^/dr = eF^v^ + / M . This is incorrect 
because it fails to obey the condition that the velocity is orthogonal to the 
acceleration.] One may correct this by hand, by generalizing / M to 

F = ^[{F aP S^r - {FapS^^v*]. (9) 
This result will be verified by starting from the action 



I = -mJd8-eJ A a dx° - -j- f d^xF^F^ - ^- J dsF^S^ (10) 

where the particle moves along the trajectory parametrized by s. This is the usual 
action of electromagnetism with the additional spin interaction added. The equation 
of motion may be obtained by considering variations with respect to the path. Using 



6S »u = s^J x ° (11) 

we have 

^ = eF>"v + (12) 

where jf is given by (Q. After the variation, we take s to be the proper time r. 
This equation may also be written as 

4- (Vm) = eF^v a + {*F af3 S a v^. (13) 
dr 

This form shows that one may consider rh = m +*F a pS a v 13 as an effective mass, 
which can be seen directly from the action. 

In the low velocity limit one may show that, calling f" — > /, 

f = V(fi-(B-vxE)). (14) 

This is a necessary result, showing that the action reduces to the correct limiting 
form. The field equations that follow from this action are discussed in the Appendix. 



1 The identity v a v a = 1 implies that v a dv a /dr = 0. 

2 An important choice has been made here. One may use F^S^ — 2*F^S fJ 'v L ' and then consider 
58^ = S^Sx 17 . In the low velocity limit this predicts that the spin force is given by / = -|(£x/j) /c, 
which is the incorrect limit. The correct limit is given by (|14p . 
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3 Generalized BMT equation 



The BMT equation starts from [12 



The conventional BMT equation is obtained by using (fT5l) with / M = 0. To obtain 
the generalized BMT equation we retain jf and obtain 

dS^ = £e_ Fa p s + 9^r _ 2 )v a S a vaF a(3 - v a Saf/m. (16) 
dr 2m p 2m vy ; p PJ ' K J 

It is helpful to write the equations in nondimensional form. We let 

x ii _^ x^/L and t — > tc/L . For example, if we consider a plane wave, 

E = E cos(kz — ut)x —> E cos(z — t)x if L = X/2tt. Also, S — > 5*//i and 

— > F^/E. Then the equations of motion for an electron, (g = 2), become, 



and 



d,v^ 

av a F^ + A/ M (17) 



dS^ 

— = aS a F^-Av< i S a r (18) 

where 

= (F af3 S a ^ - ^(F^V). (19) 

All terms are dimensionless, and the strengths are determined by the 
dimensionless constants (bringing back c), 

e FL eEh 

a = and A = — -. (20) 

mcr m z c A 

The ratio A/a = h/(mcL) shows that the spin effects on the translational 
equation of motion are generally small but, for an electron mass, become 
appreciable for hard x-rays and important for gamma rays, and may even be 
dominate for high energy gamma rays such as observed by EGRET. 
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4 Radiation Reaction 



As mentioned in the introduction, intensities of 10 22 W cm -2 have been reached, 
and this number is expected to go even higher. At such extreme conditions, the 
radiation reaction force is expected to become important, and we will examine the 
intensities where the onset of the self force becomes important before comparisons 
with the spin is made. 

It is helpful to begin with the relativistic equation of motion, ([6]), and 
consider an electromagnetic wave of the form 

E = Eh(z-t)x (21) 

and 



B = Eh(z-t)y. (22) 

This represents a plane wave of amplitude E, polarized in the x direction, described 
by any dimensionless function h. With this (jSJ) gives, 

— = ahv 1 (23) 

GST 

^ = ah(v° - v 3 ) (24) 
dr 

dv 2 , i 

= ° (25) 

^ = ahv l . (26) 
dr 

We note that ([23]) and (|2g) imply, 



v° = 1 + v 3 (27) 

which leaves the pair 

dv° 

^- = ahv 1 (28) 
dr 
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ah. 



(29) 



These imply, 



v° = 1 + (v 1 ) 2 /2. 



(30) 



Using the integral of (TJTj) in the right hand side of fl29|) yields 



v 1 = a / h{—r)dr. 



(31) 



a fascinating result. It states that the x component of the four velocity is essentially 
equal to the nonrelativistic three velocity evaluated at the proper time r. With this, 
(1271) . and (130]) . the relativistic solution is completely determined in terms of the 
nonrelativistic solution for an arbitrary ID wave form h(z — t). Having an analytical 
expression for the velocity is useful for looking at both the self force and the spin 
force in an iterative approach. 

The Lorentz Abraham Dirac equation of motion, which includes the self 
force, is[] 



where the dimensionless parameter b = ct /L, r = 2e 2 /(3mc 3 ) and the overdots 
imply differentiation with respect to the proper time. The LAD equation fell under 
bad times due to the the runaway solutions, pre-acceleration issues, or the apparent 
need to invoke non-zero size particles. [TT] For a nice introduction to the issues, with 
historical notes, one may consult Rohrlich[13], which has many of the seminal 
references and discusses the distinction between "self forces" and "radiation reaction 
forces." More recent work considers the problem in various dimensions, [14j[15j the 
effect of the magnetic dipole,[16jp2], connections to OKI ) 18 and vacuum 
polarization, [20], mass conversion |19j. and hydrogenic orbits. [21 J 

The Landau Lifschitz Rohrlich equation is obtained by starting with fl32]) but 

using 

3 We assume v a v a = 1. In the literature, some take v a v a = —1, which changes signs in the self 
force. 



aF^v^ + b{v a + v"v u v a ) 



(32) 



dr 
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dr 

in the terms multiplied by b, which leads to 



d,if 

aF CT % (33) 



dr 
where 



dif 

aF^Vp + br% (34) 



r ° = a F^ V(t + a 2 (F^F 7 % + F^v^Fjv^). (35) 



We may now compare this equation to ffTTj) . Since a » b, we may look at the ratio 
of self over spin force (ratio of constants only) 

ba 2 4 

i? = — = -aa (36) 

which shows for large fields, for any given wavelength, eventually the radiation 
reaction dominates. In fact this ratio becomes unity at / ~ 10 21 W cm -2 for 
L = 5 x 10 -5 . Figure [1] shows R on the vertical function of intensity in 

watts cm -2 and n = X/L. This also shows that, for a fixed intensity, the spin effects 
become important for short wavelengths. 

5 Solutions 

5.1 Self force 

It is helpful to begin by examining the effect of any force, JP 7 , that is multiplied by 
a small parameter A, 

dtf 

J- = a F° v v v + kT° (37) 
dr 

and to expand the solution in terms of A, 



v a = v a + A( lV a ) + ... (38) 

where ov a is the solution with T a = 0, iv a is O(A), and so on. Then, for the plane 
polarized field used above, the zero order solutions are given by (I2T1) . (|28|) . and (1291) . 
and the 0{b) equations are, 
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Figure 1: The ratio R as a function of wavelength and intensity, with n = X/L. 



d,v° 

ah( lV l )+F\w) + 0(b 2 ) (39) 

ah(( lV °) - lV 3 ) + r\ v) + 0(b 2 ) (40) 

ah( lV l )+F\ v) + O(b 2 ). (41) 
dr 

Now let us suppose AJF CT = br a , which is defined by (135]) . With this, it can be 
seen that if — i v 3 = 1 + 0(b 2 ). Calling = of 1 and dropping the subscripts so that 
iv a — > v a we have, 

dv° 

— = ahv 1 + ah(j) - a 2 h 2 (j) 2 /2 + 0{b 2 ) (42) 
dr 

and 



dr 

d\v x 
dr 

d\V 3 
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Figure 2: v z plotted against the proper time Figure 3: dz/dt vs. t, corresponding to Fi£ 
for I = 10 23 W cm" 2 . [2 



dv 3 
dr 



dv 1 

-^ = ah + ah- a 2 h 2 <f) + OVb 2 ) 
dr 



ahv 1 + ah<f) + a 2 h 2 {l - 2 /2) + Oib 2 



(43) 



(44) 



As an example, let us consider a pulse of soft x-rays, where we take 

p -((z-t)/ W f 



h 



W 



cos(£l(z — t)) 



(45) 



where the dimensionless Q determines the frequency. To zero order in b we have, 



1 

r 



1 ^ae-^ w2n2 V^h + eii(--^) + erf 



' t iwVt 



f t iwQ 



(46) 



2 J \w 2 

Before considering x-rays it is interesting to consider the case that Q = 1 = w at 
10 23 W cm -2 . This is just below the radiation reaction "threshold," but it is 
interesting to see how strongly relativistic the solution is. This is evident in Figs. |2] 
and [31 which show the four velocity as a function of proper time and the 
corresponding velocity (dz/dt) versus t (t is (roughly) laser periods, one may 
multiply by L/c to obtain time in seconds). For soft x-rays ~ 100 (50A radiation) 
and we take w = 1/10 (~ 200as pulse). The results are plotted in Figs. H] and [5] for 
I = 5 x 10 23 W cm" 2 . 

For the simple case of an infinite plane wave a typical solution is 



— asin(z — t) 



b l-asin(z - t) + acos(z - t) - — cos 3 (z-i)J +0(b 2 ), (47) 
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Figure 4: t> plotted against lOt with no radi- Figure 5: V\ 

ation reaction (solid line), and solution 
to 0{b) (dotted line). 

which will be used to compare the effects to spin. 
5.2 Spin 

Consider first the case that A/ a << 1 , valid for all but very short wavelength 
radiation. The formalism of the previous section may be used and we find in lowest 
order, 



and 



dr 



dr 



(48) 



aS„F^, 



for the same field that was used previously We already found ([271) - 
shows that S° — S 3 = const = s and 



(49) 

2ED, and gHD 



and 



S 1 = as sin(z — t) 



S 3 = ^sm 2 (z-t), 



(50) 



(51) 



and S 2 = S? , the incident value of the y component. This result may be used to 
evaluate f 1 in an iterative manner. Using the above results for the spin we find, for 
example, 



-asin(,2 — t) — A [asm(z — t) + AS^ n asin2(z — t)J . 
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(52) 



These solutions correspond to somewhat specific incident fields, but the main 
purpose is to examine the relative weight of the effects of radiation reaction spin 
forces. Comparing (J52"j) to ( pE7j) explicitly verifies the result (156]) . 

In summary, the classical equations of motion for a particle with an intrinsic 
magnetic moment in an electromagnetic field have been derived. They were derived 
from an action principle in Minkowski spacetime along with the field equations. The 
generalized BMT equation was given, and the non-dimensionalized form of the 
coupled equations was given, a simple approximate solution was exhibited. The 
effects of the spin force and the self force were considered for high intensity 
radiation, and conditions where the radiation forces become larger than spin effects 
were examined. 

The results derived here are not readily amenable to be tested in current 
laboratory experiments, but they tell us when spin is important, the onset of 
radiation reaction, and the energies and wavelengths for which self forces begin to 
dominate over spin. These results will be used to guide future work that will 
calculate observable characteristics, such as radiation emission spectra, of 
ultrarelativistic particle motion. Another interesting area with this formalism can 
be applied is the motion of neutral particles with spin in strong magnetic fields, 
such as the behavior of neutrons near magnetars. 

6 Appendix 

The action has been generalized to describe the interaction of the electromagnetic 
field with a point dipole, but we should bear in mind that the field equations also 
follow from the action. In this case we hold the trajectory fixed while considering 
variation in A^. This gives, 




(53) 



where we find the usual 




(54) 



and, using S = 5 4 



(x — x(t)) 
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i v =^-jdxp(e^SJ)^. (55) 

Thus, the equations of electrodynamics have been altered by the additional 
source i u , and it is important to consider the ramifications of any modification to 
such a successful theory as electrodynamics. First, it is easy to see that since 
F^ u = 0, {f + i>*) )lt = 0. The definitions (JUj) and ([55]) show that the four 
divergence of each of these vanish. This we have the usual conservation of charge. 

To see the effect of the new source term, consider the nonrelativistic limit, 
and consider the source defined over some small region of space. A mulitpole 
expansion of the vector potential is 

AJx) = - [ J n (x')d 3 x' + —,■ I x'J n (x')d 3 x' + ... (56) 

X J X A J 

where J n = j n + i n . Now consider some small volume where the contribution from 
j n is negligible (such as the size of a particle). Related calculations may be found 
elsewhere. |22j and it may shown that the first term gives no contribution and the 
second leads to,[23j 

* - ^ (") 

This new term actually solves a fundamental flaw in electrodynamics. The 
electron has a magnetic dipole moment which gives rise to a magnetic dipole field. 
This field can be measured classically, e. g., be detecting the magnetic field of any 
substance in which the field arises solely from the intrinsic moment (and not from 
the orbital angular momentum). As is well known, it is impossible to account for 
this observed field from j M alone: Any charge distribution confined to particle sizes 
would have to involve linear speeds far in excess of the speed of light. However, ([5 
shows that the new source term gives precisely what is observed^ 

As far as QED is concerned, we still assume that the Lagrangian density is 



C = ^ - A»d„ - mty - ~ V"". (58) 



This gives 



4 This does not predict the magnetic moment, though. That was put in by hand in (jHJ), ([5]) shows 
that g, the gyromagnetic ratio, was taken to be 2. 
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F» v v = eip^ijj (59) 

which shows that, quantum mechanically, the interaction is unaffected by the new 
spin term, although quantum spin effects might manifest themselves on other 
ways. 
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